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Abstract. This work presents a many-fermion Hamiltonian with the following 
properties: 1) is exactly solvable, 2) has a second order insulator-metal quantum 
phase transition, 3) has a well defined mean field approximation and 4) its mean- 
field ground state displays a liquid-solid transition. The phenomenon of symmetry 
breaking in fermionic self-consistent models is discussed in the light of these 
remarkable properties of the many-body model. 



Lessons from an exactly solvable Many-Body model 



2 



PACS numbers: 71.10.-w,71.10.Pm,71.10.Hf 

1. Introduction 

Spontaneous translational symmetry breaking, i.e., crystallization, is one of the most 
interesting problems in condensed matter. But as stressed by Robert Laughlin in his 
book pQ, the liquid- solid transition remains one of the outstanding unsolved problems 
in theoretical physics. To date, there is no quantum mechanical model in which 
crystallization can be observed first hand. 

In contradistinction, translational symmetry breaking occurs quite often in 
the self-consistent models of the condensed matter. For repulsive interactions, 
Ref. [2] showed that the energy levels in the unrestricted Hartree-Fock approximation 
are always fully populated. This result automatically implies symmetry breaking 
whenever the last occupied level of the translational invariant Hartree-Fock 
Hamiltonian is only partially populated. Ref. [3] proved the existence of symmetry 
breaking in the fermionic Hartree approximation for short, attractive interactions. 
The precise conditions and the mechanism of symmetry breaking in the Kohn- 
Sham equations was discussed in Ref. jl]. The Wigner crystallization of the 
electron liquid was studied in the Refs. [5] through a combination Monte-Carlo and 
Density Functional Theory calculations. The symmetry breaking in the mean field 
approximation of the bossonic Hubbard model was recently discussed in Ref. [BJ. The 
modern theory of freezing [8l |7] is based on the assumption that the linear density 
response equation of the liquid displays a crystalized self-consistent solution as soon 
as the system crosses the liquid-solid phase boundary. And the list can continue. 

The crystallization in fermionic self-consistent models, such as Hartree, Hartree- 
Fock and Kohn-Sham, is triggered by a strong coupling between the bare electrons and 
holes at the opposite sites of the Fermi surface. In well defined conditions, these models 
were shown to display robust fixed points that break the translational symmetry 
of the many-body Hamiltonian [3J H]. Almost like a rule, for metallic systems, the 
crystallization was observed to be accompanied by a gap opening in the mean-field 
spectrum. 

In this work, we deal with the translational symmetry breaking in self-interacting 
1 dimensional (ID) quantum liquids. According to a result by Mermin [9], true 
crystallization cannot occur in 1 and 2D for specific short range interacting systems. 
This result, however, does not entirely exclude crystallization in ID, as shown by 
explicit ID models [TU], or mathematically rigorous statements [TTJ 021 [T3] . ^ n an y 
case, the issue is extremely puzzling because the self-consistent models have a tendency 
to display robust crystallization precisely in lower dimensions. 

Here we propose a many-body Hamiltonian, which is a band insulator with a 
particular two-body interaction that takes into account only the electron- hole coupling 
mentioned above. While the model is fairly simple, it still displays interesting features 
such as a second order singularity of the ground state energy as function of interaction 
strength, an insulator-metal phase transition and also crystallization in its mean- 
field approximation. Based on these features, we can explicitly sec how how the 
crystallization seen in the mean-fied approximation relates to the exact many-body 
solution. 
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2. The many-body model 

We start with several considerations that will allow us to place the many-body model 
introduced here relative to the Luttingcr liquid [14 . The Luttinger liquid concept 
is now generally accepted to apply to all conducting spinless fermion systems in one 
dimension |14j . whenever such conducting behavior can be established. The Luttinger 
liquid state is, however, known to be unstable against an insulating pinned charge- 
density-wave state pT] and in fact the full phase diagram of the ID Fermi systems is 
not presently known. It will be interesting to establish a solvable many-body model 
that is relevant to a region of the phase diagram not covered by the Luttinger liquid 
concept. We will try to argue that this is the case for the present model. 

Let us start from the ID spinless Fermi gas on a circle of length L, self-interacting 
via a two-body potential v. The dynamics of the gas is generated by the following 
general Hamiltonian: 

k k 

where a k creats a fermion in the state 4> k (x) = e lkx /y/L, k = 2nir/L, n = 0, ±1, . . ., 
and p k = J2k' a k'+k ak ' n The Fermi wave-vector, kp, is assumed equal to one of the 
k n vectors. We follow Ref. [T5] and introduce the Luttinger variables as: 

ai q = a q+kp , g G [— fcj?,oo], a 2q = a q ^ kF , q e[-oo, kp). (2) 

Note that no fictitious states have been introduced. With these new variables, the 
kinetic energy becomes 

T= ^2 q(k F + q/2)a\ q a lq - ^ q(k F -q/2)al q a 2q (3) 

q> — k F q<k F 

and 



k F -k 

Pk= a lq+k a lq+ a l q +k a 2q+ a U+k-k F a 2q+k F - (4) 

q= — k F q= — co q= — k 

Our observation is that, in the limit kp—too, the third therm in Eq. Q commutes 
with all a\ q and a 2q - Mathematically, this is equivalent to say that it goes weakly to 
zero. Consequently, pk — * Pik + f>2k in the limit kp — > oo, with 

oo oo 

Plk = Y a lq+k a ^ P2k = Y a 2q+k a 2q, (5) 
q— — oo q— — oo 

being the classic density operators appearing in the Luttinger model. We can 
immediately conclude that the Luttinger model [IB] . 

H C = Y l( a lq a lq ~ a 2q a 2q) + \ Y V k{plk + P2k)(pl-k + p2-k), (6) 
q k 

can be viewed as the following limit of the ID fermi gas: 

H c = lim ^-[T+kpV]. (7) 

k F — >oo fvp 

According to this observation, the Luttinger liquid concept is for sure relevant in the 
limit of large particle densities (large kp). But in this limit, the Fourier component 
V2k F of any well behaved interaction goes to zero, and in fact the exactly solvable 
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Luttinger model excludes the coupling between the states near ±kp, contained in the 
third term of Eq. Q . This coupling is precisely at the origin of the charge-density- 
wave instability. 

Our many-body model involves the extreme case of a two-body potential, 

Vk = V [8k,2k F + 4 -2k F ), (8) 

which zooms into the coupling of the 2kp spaced one particle states and ignores 
anything else. Only the case vq<0 will be examined in this work, which is interesting 
for the symmetry breaking problem. We also restrict the one particle Hilbert states to 
fc's within the intervals [— 2kp, 2kp). Therefore, we ignore in the original Hamiltonian 
Eq. (Il| all a,k and a\ with index outside the interval [— 2kp,2kp). With the notation 
p^ =~^[pk ± p-fc], the interaction potential of Eq. [l] becomes 

v - ^[(pitir - fain o) 

Given our constraints on the k wavenumber, the operators p^fep mv °l ve only the third 
term of pk in the expansion of Eq. Q, i.e., exactly the term that is neglected in the 
exactly solvable Luttinger model. 

Before we assemble the total Hamiltonian, we need to discuss the one-particle 
Hamiltonian: 

H = ^2 ( £ k - £F)a\a k . (10) 

-2k F <k<2k F 

We would like to start from an insulating state, i.e., exactly opposite to the Luttinger 
liquid regime. We propose the following dispersion: 

f e F - A for |fc| < k F 

£ k= \ A , I, I , (11) 

[ e F + A for |fe| > k Fl 

in which case Hq describes an insulator with two non-dispersive bands, viewed in 
an extended Brillouin zone, separated by a gap 2A. The proposed full many-body 
Hamiltonian is 

H= £ (^-^)4a fe + ^[(4£) 2 -(4i) 2 ]- (12) 

~2k F <k<2k F 

To summarize, the above Hamiltonian assumes an energy cut-off at 2kp 7 or a two band 
approximation, non-dispersive bands separated by a gap 2A, and a singular two-body 
interaction that takes into account only the electron-hole couplings that were argued 
to be relevant for the translational symmetry breaking. 

3. Diagonalizing the Hamiltonian 

It is convenient to render k from ±kp and express in terms of creation and 
destruction operators with respect to the ground state of Ho, for which all the 
states with k G [— kp,kp) are occupied: 

f b q , -k F <q<0 ( c\, -k F <q<0 

a q -k F = \ , , a q +k F = { (13) 

[ 4, < q < kp { b q , 0<q<kp. 

In terms of the new creation and destruction operators, 

& = \Y.H c l+ c ^ & = l y ExM4- c * b *i> ( 14 ) 

9 9 
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where — kp < q < kp and \{q) = — 1 for — kp < q < and 1 for < q < kp. The 
interaction potential can be explicitly diagonalized by considering 

Li = P&l, I* = -ip£ p , I* = l E X{<l)[b\b q -c q c\}. (15) 

—kp<q<kF 

A simple algebra shows that [L^ Lj] = ie^L^ and that 

V= V j-[L 2 -L 2 ]. (16) 

The eigcnstates of V are the usual \ JM) states, where the highest weight states \ JJ) 
are in general not unique, and one can show that J max = \Nq, with Ao the number of 
particles. N is an even number because of our careful choice — 2kp < k < 2kp. The 
following eigenvectors can be computed explicitly: 

I Jmax, Jmax) = JJ &J4* , I Jmax, -Jmax) = JJ b\ C \^ ,{^) 

0<q<k F -k F <q<0 

and the manifold | J max M) can be shown to be non-degenerate (i.e. there is just one 
highest weight vector | JmaxJmax)- The lowest energy state of V is | J m axO), which can 
be obtained from the above vectors by applying the Lzp operators ^Nq times. 

We now consider the whole Hamiltonian H=Hq+V . The new representation of 
Ho is 

H ^Aj2[blb q -c q c g }. (18) 

H commutes with L 3 but not with L 2 . It is interesting to notice that there is a 
competition between V, which favors electron-hole pair formations, and Ho, which 
does not. As we shall see, this competition will ultimately lead to a quantum phase 
transition. The model can be solved by employing the SO(4) Lie algebra. Indeed, if 
we define 

Ki = \Y,- kF < q<kF X{q)[b\c\ + c q b q ] 

K2 = -hT- kF < q<kF W\-c q b q ] (19) 

K 3 = \H- kF < q<kF \ b \ b q - CqC\], 

then 

[Li,Lj] = [Ki, Kj] = ie ljk L k , [Li,Kj] = ie ijk K kl (20) 
precisely the SO(4) Lie algebra. A straightforward calculation gives 

H = 2AK 3 + V ^[L 2 -L 2 }. (21) 

If A ee \{L + K) and Y = \{L - K), then 

[X i ,X j ]=ie ijk X k , [Y i ,Y j ]=ie ijk Y k , [X i ,Y j ]=Q (22) 

and 

[H, A 2 ] = [H, Y 2 } = [H, A 3 + Y 3 ] = 0. (23) 

The operators A and Y commute with the particle number operator, therefore we can 
restrict the discussion to the quantum states with precisely N particles. 

Let us denote the common eigenstates of A 2 and A 3 by \jm) and the common 
eigenstates of Y 2 and Y 3 by \j'm'). We verified that the maximum values of j and 



Lessons from an exactly solvable Many-Body model 



f are jmax=imax = i-^o- Now, for fixed j, j' and M, the vector space spanned by 
<jMi M ^ — \jm) ® \ j'M — to) with to taking all allowed values, is invariant for H: 

H(p m = y(a m _i^„-i - 2a m (/) m + a m+ i(/> m +i) + v m 4> m , (24) 

L 2 2 

where, we dropped the upper indices to ease the notation. The coefficients appearing 
above are given by: 



a m = + 1/2)2 _ TO 2][ (/ + 1/2)2 _ ( m _ M) 2] 

v m = f + 1) + f(f + 1) - M 2 + 2m(M - to)] (25) 
+2A(2to — M) + 2^a m . 
The vector 4> m must be set to zero if |m| > j or \M — m\ > j' ', a statement that 



clarifies the allowed values of to. Eq. (24 1 allows one to calculated the whole energy 



spectrum of H and solving it is no more complicated than diaganolizing a one-particle 
tight-binding model in ID. The many-body eigenvalues fall into distinct manifolds 
that can be labeled by j, j' and M. To calculate thermodynamic functions, we also 
need to calculate the degeneracy of each manifold and for that we need to compute 
how many highest weight vectors \jj) <g> \j'j f ) are there for each j and j' . This will 
not be done here. 

4. Thermodynamic limit 

Things greatly simplify in the thermodynamic limit L— >oo, when j, j' and M take 
macroscopic values (proportional to L). For large L, we normalize to by j and work 
with x = m/j (\x\ < 1) as our variable, which now becomes continuous. For fixed j, 
j' and M and with the representation ^ = c m <j) m , where c becomes a function of 
x, the action of the Hamiltonian per unit length becomes: 

—Hc(x) = -^rd x a(x)d x c(x) + v(x)c(x). (26) 
L L l 

For each manifold {j,j',M}, the potentials a(x) and v(x) can be easily derived from 



Eq. (251 



v(x) = v {j/Lf [1 + (j'/j) 2 - (x - M/j) 2 - x 2 ] 

+ 2A(j/L)[2x - M/j] + 2v (j/L) 2 a(x). (27) 
Eq. (26 1 is defined on the interval where the factor under the square root in a(x) is 



positive and zero boundary must be imposed at the ends of this interval. 



From Eq. (26 1, one can see that, when L—>oo, the energy per unit length is given 
by v(x). The ground state is contained in the manifold j = j' = j max = 7N0, and 
M = 0. This can be established analytically but we have also verified the statement 
numerically. For this manifold, v(x) is equal to: 

v„ lin (x)=nlv (l-x 2 )/4: + n Ax, (28) 

with n = lim N a /L — kp/n. Given that x is constrained to \x\ < 1, the minimum 

L — >oo 

of this potential, which defines the ground state energy per unit length, corresponds 
to xq — max{ 2 ^ g , — f } and the energy per unit length is: 

E n { -«oA, 2A > n \v \ 
lim = <^ . (29) 

L \ "oA[^ + ^] 5 2A<no|«o|. 
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Thus, in the thermodynamic limit, there is a sharp change in the behavior of the 
ground state energy as a function of the model's paramaters. We have verified 



numerically that, indeed, Eq/L converges to the values predicted in Eq. 29 

To get more insight, we calculated the expectation value of Hq (= 2AK3) on the 
ground state of H, divided by L: 

(H ) f -«o A, 2A > n \v \ 
Urn = { (30 
L \ 2A>o, 2A<n K|. 

From here, we conclude that for vq larger than the critical value v c — 2A/ti , we have 
macroscopic occupation of states with \k\ > hp- 

5. Mean field analysis 

We define the mean field problem by the following substitution: 

v^^K&piZ -&& ( 31 ) 

where (p^ ) denotes the expectation value of p^) on the ground state of the mean 
field Hamiltonian, which has to be calculated self-consistently. If a=(p2k F ), then 
(^2fe )=Re[a] and (p\^ F )=ihn{a\. We assume in the following that a is a real positive 
number. The mean field Hamiltonian becomes: 

H MF = J2 {Mb% - c,ct) + ^°(6t c t + Cqhq) } . (32 ) 
9 

We can diagonalize the quadratic mean field Hamiltonian by using the following 
Bogoliubov substitution: 

b q = cos 8bq + sin 9c q: c q = — sin 9b q + cos 9c qi (33) 

with 

tan20 = ~^, 0e[0,n/2]. (34) 

The mean field Hamiltonian becomes 

A 
cos26> 

9 

and 

1 sin 28 x—i _ cos 29 



^ = -^Efe+4 5 «- 1 )' (35) 



P2k f 



E(^-c g ct) + ^^E(6tgt + c 9 6 9 ). (36) 



2 Z^v-9 9 -v-q/ 1 2 
9 



It remains to determine a = (p2k F )> self-consitently. Of course, there is always the 
trivial solution a—0 but we will show that, if vq is negative enough, the mean field 
approximation has nontrivial self-consistent solutions. 

Indeed, let us assume 1 > cos 29 > 0, in which case the ground state of H MF has 
all the b q and c q states empty and the expectation value of p^ k becomes 

a = VaiVoshi20, (37) 



which together with Eq. (34 1 leads to 



cos20 = -— anda= ^/l-f— V- (38) 
n v 2 V \n v J 
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Figure 1. (a) The energy gap of the many-body Hamiltonian H as function of 
no vq. (b) The energy gap between the first and second levels of the manifold 
{imaxJmaxAf = 0} (solid line) and the energy gap of the mean field Hamiltonian 
-Hmf (dashed line). Both are plotted as functions of nrji>o- The only input for 
these calculations was iVo=4000. 



This self-consistent solution is in line with the starting assumptions that a > and 
cos 26 > 0, therefore it is a valid solution. This solution exists as long as 2A < no|fo|, 
otherwise the the cosine will exceed its maximum value of 1. 

The mean field approximation predicts a ground state energy per unit length 

E mf ( -n A, 2A > n |u | 
Urn -5- = { (39) 
L ^°° L { nlvo/2, 2A < no|«o|, 

which also displays a sharp transition. 



6. Discussion 



The first observation is that both the exact and the mean field treatments predict the 
same critical value of v , where the character of the ground state changes. But there 
are several differences between the two. The exact ground state energy per unit length 
Eq/L and its first derivative with respect to Vq are smooth as v crosses the critical 
value. The second derivative, however, has a jump. For the mean field approximation, 
already the first derivative of E^ F / L with respect to vo has a jump at v c = — 2A/rto- 

But the most interesting difference concerns the gap of the energy spectrum. In 
Fig. [lja) we plot the gap of H (no length normalization) as function of uqVq, with 
No fixed at 4000. The gap of H is given by the difference between the ground state 
energies of the manifolds {j max , jmax, M = ±1} and {j max , jmax, M = 0}. As one 
can see, the exact solution indicates a transition to an un-gapped system. In fact, 
above v c , the ground states of the manifolds {jmax, Jmax, \M\ > 0} form a continuous 
sequence of energy levels. Therefore, at v c we have a true insulator-metal transition. 

As already mentioned in the Introduction, the non-trivial mean field solutions 
found so far predict metal-insulator transitions at the symmetry breaking. The 
situation is similar in the present case, even though we start from a band insulator. 
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In Fig. [TJb) we plot the gap of H MF as function of uqVq. We see the mean field gap 
first decreasing, but never reaching zero, and then increasing robustly. Therefore, the 
mean field ground state becomes more insulating as uqVq is increased. At first sight, 
this seems to be in total contradiction with the insulator-metal transition of the exact 
ground state. But if we look at the gap within one manifold, we actually see a very 
strong correlation. In Fig. [ljb) we also plot the difference between the second and 
first eigenvalues within the {j m ax,imax,JW = 0} manifold. As one can see, the mean 
field gap actually reproduces the manifold gap quite precisely as UqVq gets larger. 

The non-trivial self-consistent solution of the mean field problem implies 
crystallization. Indeed, the nontrivial expectation value of (0 M p|p(^)|0mf) = 
j J2k e~ lka '(pk) mf leads to a modulated particle density: 

/ 2 A v 2 1 



nl.r] = />„ j l + cos(2fc F a;)yi- ^— — j j . (4()i 

It is easy to check that if a is a self-consistent solution, then e lip a is a self-consistent 
solution too, and the only change in n(x) is the appearance of the phase factor 
inside the cosine: cos(2kpx)^cos(2kFX + <p). The expectation value of p(x) on the 
exact ground state, however, is independent of x and equal to no. In other words, 
crystallization is absent in the exact ground state. But as already pointed out, the 
system becomes gapless so it will be interesting to see if crystallization can occur in 
the low lying excited states. 

To answer the last question, we start from the action of p^. '■ 

P2k F ( Pm — A M,m C Pm±l ' l 4i J 

with A^ M = y/(j Tm)(j ±m + l) + =F (M - m))(j' ± (M - m) + 1). In other 

words, takes the {if M} manifold into the {jj'M + 2} manifold and p^' p takes 
the {jj'M} manifold into the {jj'M — 2} manifold. Given these actions, and the 
fact that the ground state energies of the manifolds {jmaxjmax^ = ±1} are the same 
and converging to the absolute ground state energy in the thermodynamic limit, we 
reached the conclusion that the mean field ground state actually relates to the states: 

=" X) + ^- j - M=+1 \ (42) 

where M denotes the ground state within the manifold {jj'M}. 

In Fig. [2] we present plots of the eigenvectors discussed above. Although Eq. [26] 
suggests that H becomes just a multiplicative potential in the thermodynamics, the 
eigenvectors actually display a quantum spread even for large L. This is because the 
potential v m becomes more and more flat as L is taken to infinity. In Fig. [2ja) we 
show the potential v m for the manifolds {jmaxjmax-^ = 0,±1}. The inset shows a 
global picture of v m for uqVq — — 5A and manifold {j ma , x j m axM = 0}. The main figure 
is a zoom into the bottom of the potential. One should notice the very small scale on 
the vertical axis, which shows how flat the potentials are. One should also notice that 
the minimum of the potentials occurs at certain values of to, that shift to right as M 
is increased. Although the potential of the manifold {jmaxjmax-^ = 0} is lower, this 
cannot be distinguished in the figure because the differences are miniscule. 

In Fig.^b) we show the coefficients c, n of the expansion M ^=J2 m Cmfim M ' '■ 
Again, the inset gives a global picture of the eigenvectors, while the main figure zooms 
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Figure 2. (a) The potential v m corresponding to the manifolds j = j' = j ma x and 
M = 0, ±1. The inset gives a global view of v m for M = 0. (b) The amplitudes c m 
of the ground states >&o corresponding to the manifolds {(jmaxjmaxM = 0, ±1)}. 
The inset gives a global view of f° r {(imaximaxM = 0}. The input parameters 
were 7Vq=4000 and novo = — 5A. 



into the region where c m 's take large values. One should notice the considerable 
spread of the eigenvectors and that the maximum value of c m occurs at certain 
m's that drift to the right by about one unit when M is increased by two units. 
This drift is essential in order to see large expectation values of p^ , since these 

operators increase/decrease m by exactly one unit when applied on 4>m • And 

indeed, we have verified numerically that the expectation value of P2k F —P2k F +P2k F 
on the vector written in Eq. [42] is consistent with the value of a provided by the mean 
field approximation, within 3 significant decimals. 

7. Conclusions 

The exactly solvable ID many-fermion model have given us unprecedented insight 
into the phenomenon of translational symmetry breaking observed in mean field 
approximations. The mean field of the present model has a robust solution displaying 
the crystallization of the ID fermions, provided the strength of the interaction exceeds 
a certain threshold. The crystallization is absent in the exact ground state, but can 
be found in the first excited doubly degenerate level, whose energy converges to the 
ground state energy in the thermodynamic limit. Within this first excited energy 
level, we have constructed, explicitly, a many-body vector which displays exactly the 
same particle density as the mean field ground state, therefore establishing a precise 
connection between the mean field and exact solutions. 

We have also demonstrated that the model exhibits a band insulator-metal 
transition. An intriguing question arises. Are we witnessing the transition from a 
band insulator into the Luttingcr liquid? According to the general accepted view, 
the Luttinger liquid state is relevant to all un-gapped ID fermions, so one would be 
inclined to answer positively to this question. But the model presented here is in many 
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respects complementary to the Luttingcr model, therefore we are inclined to believe 
that the metallic state is rather coupled to the charge-density-wave state. This issue 
will be investigated in the near future. 

At the end, we should mention that there are hopes that the exactly solvable 
model is stable enough to allow perturbative extensions that touch on realistic systems. 
For example, replacing the two perfectly flat bands with weakly dispersive bands 
should not pose major problems. Small variations of Vk seem also to pose no major 
difficulties, but it is not clear at this point if such variations can be pushed to the limit 
of realistic interactions. Clarifying these points will be important for understanding 
the impact of the model on the phase diagram of the ID fermionic systems. 
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